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The point symmetry group is studied for the generalized Webster-type equation describing non- 
linear acoustic waves in lossy channels with variable cross sections. It is shown that, for certain 
types of cross section profiles, the admitted symmetry group is extended and the invariant solutions 
^ \ corresponding to these profiles are obtained. Approximate analytic solutions to the generalized 

Webster equation are derived for channels with smoothly varying cross sections and arbitrary initial 
j /^ ■ conditions. 

<n: 

I. INTRODUCTION 

Q-f The Webster equation [lH3j describes waves propagating in pipes, horns, concentrators, and other waveguides 
^ ch characterized by a varying cross section S(x): 

1 o ( d P \ 1 d 2 P 

S(x) — I — — —-r = , (1) 



S(x) dx \ dx J c 2 dt 2 

\ where x is the coordinate measured along the waveguide axis. Equation ([I]) is applicable to pipes whose characteristic 
width is small compared to the wavelength. In addition, the cross section is assumed to vary slowly along the rr-axis: 
the area S(x) varies only slightly as x varies by a quantity on the order of the pipe width Q. Note that, in the general 

1 case, Eq. ([!} cannot describe a wave propagating in the pipe in one direction (e.g., in the positive direction of the x 
C*") 1 axis): 

t-h ; 

Od' P(M) =A(x)$(t-<p(x)) , (2) 

o: 

^SJ . where $ is an arbitrary function. Indeed, substituting Eq. ([2]) in Eq. (Q}, we obtain the relations 
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da; I c 2 ' dx K ' dx V da; 



For the boundary conditions A(x = 0) — po, S(x — 0) = 1, we obtain from Eqs. ([3]) 

A=-^—, S=(l±x/R) 2 . (4) 
1 ± x/R 

Here, the plus sign corresponds to a spherically expanding pipe and, hence, a spherically divergent wave. If, in Eqs. 
(U]), we take a minus sign, we obtain a tapered pipe and a wave converging to the center of the sphere x — R. As 
R — » 00, the pipe becomes homogeneous (S(x) = const) and the wave becomes plane. 

Thus, it is only in the case of a plane or spherical wave that ([l} describes a wave traveling in one direction with 
an arbitrary profile In all other cases, the inhomogeneity of the cross section S(x) causes multiple reflections 

of signals propagating in opposite directions, which leads to the formation of a complicated wave field pattern as a 
combination of standing and traveling waves. Study of the influence of inhomogeneity S(x) on the field structure gives 
a chance to solve the inverse problem, i.e. it helps to reconstruct the unknown function S(x) using characteristics of 
the reflected (transmitted) wave. This problem was discussed in Q for small-amplitude acoustic waves in a tapered 
waveguide and also in Q for finite-amplitude waves. Solution to the inverse problem can serve as a classical illustration 
of the effect of the dimensional reduction that was discussed recently in [6|-|8( . 

However, in the case of slowly varying S(x), within distances on the order of several wavelengths, it is possible to 
neglect multiple reflections and consider a wave traveling in one direction. This statement holds for both conventional 
Webster equation (1) and its generalization to the case of a nonlinear dissipative medium filling the pipe. The 
generalized Webster-type equation appears in problems on propagation of intense sound [9|, Il0| . It differs from Eq. 
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(fTj) in that it contains two additional terms describing nonlinear and dissipative effects. We represent this equation 
in the form 

d 2 p 2 d 2 P 2 d^ n S( x ) dp _|_ £ d 2 p 2 d 3 p 
dt 2 dx 2 dx dx c 2 p dt 2 p dxdt 2 

Here, e and b are the nonlinearity and dissipation parameters (the notations are the same as those given in 

EH) 

and 

p is the density of the medium. Both the initial (U} and generalized ([5]) Webster equations are suitable for describing 
traveling and standing waves. 

In the situation where each of the terms appearing on the right-hand side of the equation is small compared to the 
terms appearing on the left-hand side, a traveling wave can be considered. In this case, using the method of a sl owly 
varying profile it is possible to reduce the order of nonlinear equation ([5]). Following the standard procedure 
we change from the variables x and t involved in Eq. ((5]) to new independent variables: the slow coordinate x\ = 5x 
(where <5 is the small parameter of the problem) and the time r = t — x/c in the coordinate system traveling with the 
velocity of sound. Ignoring small terms on the order of S n , where n ^ 2, we arrive at the evolution equation 

dp e dp b d 2 p p d 

Tx - —p P d-r ~ 2^W + 2 Tx (ln S{X)) = ° ■ (6) 

This equation is used not only as a model of wave propagation in a pipe but also for calculating the acoustic field in 
an inhomogeneous medium in the geometric acoustics approximation [a,[l3j]; in the latter case, it plays the role of the 
transfer equation represented in terms of the ray coordinates. The axis of a ray tube is the geometric ray calculated 
from the eikonal equation, and the function S(x) represents the cross section of the ray tube. 

The study of even the simplest linear Webster equation (TTJ reveals rather interesting effects, such as, e.g., tunneling 
of sound waves in the case of sound propagation through a tapered waveguide The change to generalized Web- 
ster equation (J6j> (GWE), which contains additional contributions due to nonlinear effects and absorption, makes it 
possible to study problems on the propagation of finite-amplitude sound waves in absorbing media and, in particular, 
acoustic monitoring of propagation media [5]. For these purposes, exact and approximate analytic solutions to the 
finite-amplitude acoustic wave equation acquire special importance, as in the case of the Burgers equation and its 
generalizations for homogeneous media (see [lfj and Sect. 2, Ch. 7 in [11]). 

Although solutions to the GWE simultaneously allowing for the effects of nonlinearity, absorption, and inhomogene- 
ity are important for describing the behavior of sound waves in nonlinear absorbing media, their analytic derivation 
is a difficult problem even in the case of using approximate methods. The standard practice is either to neglect 
dissipation, which makes it possible, by changing the independent variable, to reduce the initial equation to the Hopf 
equation (10|, or to assume that the nonlinearity is small and to solve the sound wave equation by the method of 
successive approximations; the latter approach has been used to analyze the second harmonic behavior in a sound 
channel with a variable cross section On the other hand, it is evident that, as the cross section of the channel 
decreases, the second harmonic amplitude grows faster than the fundamental harmonic amplitude, which necessitates 
analyzing the GWE for a finite-amplitude sound wave. This paper is devoted to finding an analytic solution to the 
GWE under these conditions. 

The paper contains five sections. The second section formulates the initial equations for the theoretical analysis of 
sound wave propagation in a medium with allowance for nonlinearity and absorption. In addition to the initial GWE, 
we consider its modified analog, which simplifies the analytic investigation by group theoretical methods. The latter 
prove to be the most effective instrument for constructing analogs of the exact analytic solutions obtained earlier 
for a homogeneous medium. In the third section, we determine the point transformation group for the modified 
GWE and show that, for specific cross section profiles, this symmetry group can be extended. We determine the 
invariant solutions corresponding to the aforementioned symmetries and compare these solutions with those obtained 
for a waveguide with a constant cross section. The fourth section is devoted to constructing an approximate analytic 
solution to the GWE for arbitrary varying cross sections and arbitrary initial conditions. The consideration is based 
on the theory of approximate transformation groups, which allows determination of an approximate symmetry and 
the corresponding approximate analytic solution. The small parameter used in the calculations is the slowness of the 
variation in the waveguide cross section profile S(x). In the fifth section, we briefly formulate the main results of our 
study. 
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II. BASIC EQUATIONS 

To analyze the nonlinear effects that accompany sound wave propagation in an absorbing medium in a channel 
with a variable cross section S(x), we use GWE (j6|), which can be represented in the form 

S-^-^ + |^( ln ^)) = °> K0,r)=P(.). (7) 

Unlike the equations discussed above, Eq. ([7]) is expressed in dimensionless form. The change from the physical 
variables involved in Eq. (j6|) to the more convenient normalized variables appearing in Eq. ((7|) is performed through 
the following substitution in Eq. ((6|): 

c T 
x-> — x , t y , pop . 

uj uj 

Here, the normalizing constants uj andpo have the meaning of the characteristic frequency and signal amplitude values, 
respectively. The two parameters involved in Eq. ([7]) are determined by the following dimensionless combinations of 
constants: 

epo buj 

a = 



c 2 p ' 2c 2 p 

Their ratio a/v is called the acoustic Reynolds number [TTj . It characterizes the relative contributions of nonlinear 
and dissipative effects to the distortion of the wave profile. When ajv is large, nonlinear effects predominate; when 
this quantity is small, dissipative effects are dominant. Without loss of generality, in Eq. Q we set S(Q) = 1. 

We can eliminate the last term from Eq. ([Jj by changing the variable x and introducing the absorption as a function 
/j, of the coordinate along the channel: 

C=/dx/v^), P VS = u, n = u^/S(^0) > (8) 



Then, in terms of the new variables, Eq. ([7]) takes the form 

du du d 2 u „ . . , . . . 

dC- au fr-^ =0 > «(0,t) = P(t). (9) 

Introducing the new variable q related to u by the formula u = 2(dq/dr), we replace Eq. (J7J) by the modified GWE 
(MGWE) 

dq /<9g\ 2 d 2 q 

The change to the new variable q in Eq. ([9]) increases the order of this equation. However, its single integration with 
respect to r yields evolutional equation (|10[) for q. This procedure determines q accurate to the function C(£) (in 
Eq. (fT0|) . this function is omitted), whose choice is fairly arbitrary. For example, for solutions to Eq. (flO|) that are 
periodic in r, this function can be chosen so as to make the period average value of q zero for any values of £. At the 
same time, it is evident that the choice of C(£) does not affect the physical meaning of u. 



^-M^z) -M^4=0' 9(P,t) = W(t). (10) 



III. SYMMETRY GROUP AND INVARIANT SOLUTIONS TO THE GENERALIZED WEBSTER 

EQUATION 

For Eq. (1101) in the case of an arbitrary inhomogeneity profile /i(£), the admitted point transformation group is 
given by three infinitesimal operators: 

d d d t d 

The first two of them represent the translation operators with respect to the variables r and q, which are evident 
from the physical point of view; the third operator corresponds to the Galilean transformation group. Transformation 
group (|11|) can be extended for the cross section profiles of a specific type: 

M ln(M(c)) ) =&(c) ' b{c) = Po + AC + ^ ' M = const ^°- ( l2 ) 
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Condition (|12p plays the role of the classifying relation separating the specific profile types for which transformation 
group (jlip is extended by the additional operator X^: 



X * = b dC + 2{ M+ T C )d-r 



Mq- 



(r 2 + 2/dC^) d 2 6\ d 



dC 2 J dq 



(13) 



Classifying relation (jT2"j) is a first-order differential equation for the function /x(C); being explicitly integrated, this 
equation determines a three-parameter (the parameters are fa, fa, and ^2) family of curves in the {Cm} space: 



In( M /i/)=d(C) = M | 



(14) 



where the form of the function d(£) depends on the relative values of the parameters fa: 



d(C) 



2M 



M i 



arctan 



/3i+2feC 



arctan 



0i 



M i (V^? -2feC)(-y/^; -4^0^+^) 

I V -3? -4/30/32 11 (V^f-4/3ofe+/3i+2,3 2 C)(V/5i 2 -4/3o/32-/3i) 



0? < 4/3 ^2 ■ 

/3 2 = 4/3 /3 2 ■ 

fa 2 > Ifafa . 



(15) 



The choice of M = corresponds to a channel with a constant cross section, d//(£)/d£ = 0; then, classifying relation 
([i"2j) is automatically satisfied for any /J^. In this case, the following three operators appear instead of the operator 
X A : 



d_ 



x 



42 



T_d_ 

2d^ 



X. 



43 



d_ 



d_ _ (r 2 + 2C,v) 9 
9t 4a <9g 



(16) 



The first of them, X41, is the translation operator along the £ axis; the second operator, X42, represents the dilation 
transformation; and X43 corresponds to the projective transformation group. In addition to operators (|16[) , for a 
channel with a constant cross section, \i = v, the MGWE also allows the infinite subgroup operator 



aq 







Ok 



d 2 k 



Xoo = k((, t) exp — , — - [i—r = 



(17) 



Here, the linear parabolic equation, which is satisfied by the function of two variables fc(£, r), can be represented in 
terms of the variables {x, r}: 

dk d 2 k 
dx <9t 2 

The latter fact will be used by us in constructing the approximate point symmetry for the MGWE. We note that the 
symmetry group given by Eqs. dill) . (I16p . and (1171) is well known in the theory of the modified Burgers equation [l5j], 
to which the MGWE is reduced in the case under consideration. 

In constructing the invariant solutions to the MGWE, we concentrate on studying solutions that are invariant 
with respect to the one-parameter group with the operator X4, because this operator represents an analog of linear 
combinations of operators (fl6)) . the use of which for channels with constant cross sections yields the well-known and 
physically meaningful particular solutions to the nonlinear Burgers equation. 

Taking into account the initial conditions for the MGWE (|10[) . we represent the desired solutions by using two 
invariants J\ and J2 and the operator X4, where 



^e- 
2b 



as Ji = W(J2) or, in an explicit form, as 



C 



-d(C') 



dC"e 



"od(C") 



■h 



:LP-«»/a 



Vb 



(18) 



q =C d ^ 



2b 



d(C) 



HO 



-rf(C') 



» P d(e") 



dC"e' 



A = 



re 



-d/2 



Vb 



(19) 
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To solve the MGWE, we substitute the latter expression in the initial equation, then we obtain a second-order ordinary 
differential equation for the function W(X): 



<w + °(drj + "" + l »-2{ir)- MW + -^ x -' > ' < 20) 

As examples, we consider several specific cases of the general relations. 
Example 1: 02 = 0, /3o,/?i 7^ 0. 

This type of a channel with a constant cross section corresponds to a combination of the translation operator along 
the x axis and the dilation group. For 02 = 0, from Eqs. (fl5|) and (|20| it follows that the solution is self-similar 
(compare with [13]); depending on the relative values of the parameters 0\ and M, it is realized for the following 
functions S(x): 

s(x ,Ul + (M + (3 1 )x/M 2M/iPl+M) , fh/Mjk-l, 
\exu(2Mx/0 o ), Pi/M = -1 . 

For an exponential dependence of the cross-sectional area on x, i.e., for 0i/M = —1, the variable A = r/V/3o does not 
depend on (; then, the solution to Eq. (f2"U|) is expressed in integral form 

\ = J dW(C exp(-2aW) + (M/2a 2 )(2aW-l)y 1/2 + d, C ,Ci=const, (22) 

and has bounded periodic solutions for M < and Co < 0. 

For a power-law dependence of the cross-sectional area on x, the self-similar variable A and the function d appearing 
in Eq. (fT9|) have the form 

r / 8i + M \ WT) M ( Si \ 

"TK^-^V • d= wA 1+ i c )- ft+M *°- (23) 

Example 2: 0i = 0, Q ,0 2 ^ 0. 

This type of a channel with a constant cross section corresponds to a combination of translation operators along the 
x axis and the projective transformation group. In our case of a channel with a varying cross section, the functions 
A and d appearing in Eq. (|19[) take the form 

Te~ d/2 M ( Hh \ 

A= n „ , d=-=arctan L/^C , /3i = 0, (24) 

where the function £ is related to cc by the formulas that implicitly determine the dependence of the cross-sectional 
area S on x: 

^tan m S ) , * = A / ^ cos - In s) . (25) 

1 2M J ' 2M J y/s \ 2M J V ' 

1 

We note that the self-similar solutions obtained for the exponential and power-law cross section variations (|2"Tj) had 
been discussed in the literature [l(|; as for integral representation (|22p of the solution to the MGWE, it has presumably 
been obtained for the first time in this paper. The invariant solution given by Eqs. (|T5)) . and (1231) has never 

been reported in the literature. 



IV. APPROXIMATE SYMMETRY GROUP AND APPROXIMATE GROUP INVARIANT SOLUTIONS 

TO THE MGWE 

The invariant solutions to the MGWE obtained in the previous section have the following drawback: being exact, 
they can only be obtained for certain specific channel profiles and initial conditions, which are given by classifying 
relation (|12[) and the solution to equation (|20p. In this section, we present alternative solutions, namely, approximate 
analytic solutions to nonlinear boundary- value problem (|10p . and these solutions can be constructed for arbitrary initial 
conditions. An instrument for constructing such solutions is the approximate symmetry group, and the condition for 
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the existence of the latter is the presence of the small parameter related to the relative slowness of variation in the 
cross-sectional area along the waveguide axis; i.e., the smallness of the derivative d(ln/x(£))/d£ = hq/ijl <C 1 (here, 
the subscript denotes the derivative with respect to the corresponding argument). In this case, the symmetry of the 
boundary-value problem under study is represented by a series expansion in powers of the small parameter, which 
allows us to obtain approximate analytic solutions to the problem with an acceptable accuracy. 

To construct an approximate analytic solution for a channel with a slowly varying cross section, we use the 
renormalization- group symmetry algorithm [16j for boundary- value problem (|10[) . According to perturbation the- 
ory, this algorithm allows us to extend the solutions in nonlinearity parameter a to the region of finite values of 
this parameter. The general description of the algorithm can be found in (l6| . and the details of the corresponding 
calculations performed in the case of solving the boundary- value problem for the modified Burgers equation are given 

in [nf. 

To take into account the transformation of parameter a, we include the latter in the list of independent variables 
and represent the infinitesimal operator of this transformation as 

*. 11:1 • (26) 



da a t 

The desired renormalization-group symmetry operator is obtained as a linear combination of operator (|26p with 
£(a) = 1 and infinite subgroup operator (|17[) . which (as we have shown in the previous section), in the zero order in 
/z^/jLt, is allowed by boundary- value problem (QI 



R=l + (k^(Cr,a)e,J-^)- q -)l. (27) 



da \ \ fi J a 

Here, the function of three variables (£, r, a) obeys linear parabolic equation (fTTj) with the initial condition 
k^(0, r, a) — W(r)/a determined by invariance of the solution at a — > with respect to the renormalization-group 
symmetry operator (|27p . As a result, we arrive at the expression 

oo 

= -K a , K{a,x,r) = [ d£e^G(z,r-0, G(x,t) = -j^e'^ . (28) 
a J y Attvx 

— oo 

Here, the subscript marking function K denotes the partial derivative with respect to the corresponding argument, 
K a = dK/da. 

Finite transformations of the continuous group are related to the infinitesimal transformation in a one-to-one manner 
by the Lie equations, i.e., equations of the characteristics for the first-order partial differential equation conjugate to 
the operator (f2~T|) . (|28p . The solution of the Lie equations for operator (|2~T1) . ([25)1 yields the following approximate 
analytic solution to initial problem (JTT 



(29) 



which is valid in a medium with a slowly varying cross section, <C 1. In fact, the derivation of solution ()29|) from 
the solution obtained for a channel with a constant cross section consists in the presence of a factor v j \i other than 
unity. 

The advantage of the renormalization group method is the possibility of sequentially refining the resulting analytic 
approximations. As applied to the problem under study, such a refinement (in the next, i.e., first-order, approximation 
in fJ,{/fJ>) is achieved as follows: the function k^((,r, a) in generator (|27]l is replaced by fcW ((, r, a), for which, instead 
of using the solution to the parabolic equation (fTT)) . we use the solution to the inhomogeneous parabolic equation 

dk& d 2 kW 

Here right-hand side is proportional to the gradient of the channel cross section fi^/fi <ti 1 and linearly depends on 
the function of the zero-order approximation in this gradient: 

A = ak^q (0) ^/^ = (/*c//*)("/a)*« ln [1 + iy/n){K - 1)] . (31) 

Equation ()30|) is obtained at the stage of calculating the renormalization-group symmetry operator (|27|) from the so- 
called group determining equation, where the contributions proportional to (which were omitted at the previous 
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step) are calculated using the zero approximation results (j2"8")l and (|2T)]) . The solution to Eq. ([3U)l yields a modified 
(due to the contribution with the cross section gradient) expression for the function k^°\a,x, r) fc^^a, X, r) 



a; oo 

fc« = - a K a - J dx' J dr'G( 



X — X , T — T 



(32) 



fi' = fi(x'), K' = K(a, x', t') 



The substitution of k^ 1 instead of in infinitesimal operator (j2"T)) and the subsequent solution of the Lie equations 
leads to a refined approximation for the desired solution: 



l( i) _ M 



(X OO 
1 + - (K - 1) - - / dx'^L f dT ' G f x — x\t — t') 
V M J M J 



1 - X' 



For small values of the nonlinearity parameter a, the first two terms of the expansion of solution 
the nonlinearity parameter have the form 



(33) 



in powers of 



<K<?> 



va 

T 



dx' 



vfj, x 



dT'G(x- X ',T-T')((K')W) 2 



Q(a 2 



(34) 



where and Ka represent the values of the partial derivatives of function K calculated in the limit a a 



substituting the periodic initial condition W(t;) 



0. By 

cos£ in K and calculating the resulting integrals, it is possible to 



verify that the expression for q pt agrees well with the result obtained in Q for a harmonic initial perturbation. 

In closing this section, we present the form of the solution to problem (fTUl) for a periodic initial condition W(^) 
cos £: it is given by Eq. (f3"3"]) with function K determined by the expression 



K = Io(a/v) + 2 Ik{a/v) cos fcr e" 



vk x 



(35) 



fc=i 



The use of Eqs. (|3"3"|) and (|3"5"|) allows us to calculate the nonlinear distortion of the spectrum of an acoustic wave 
propagating in a waveguide with a varying cross section. This opens up better possibilities for diagnostics of acoustic 
propagation paths, as compared to the weakly nonlinear limit [5j. As an example, in Fig. [T] we represent the solutions 
to the MGWE, q = versus r, the solutions being obtained from Eqs. (f3"3"]) and ([33]) for different values of the x 
coordinate along the axis of a channel with an exponentially varying cross-sectional area, [ijv = exp(air). 

To estimate the accuracy of the approximate analytic solutions obtained by us, we numerically solved initial equation 
(ITU|) . Comparison of the curves plotted with the use of analytic results (13"3"|) and (1331) with the curves obtained from 
the numerical solution of initial equation (|10l) shows good agreement between the numerical and analytic results for 
the case of a moderate nonlinearity a/v = 1 with an accuracy of up to fractions of percent. The results of calculating 
q = qW for a greater nonlinearity ajv = 10 show a difference between the results of numerical and analytic calculations; 
the difference increases with increasing distance along the pipe, as one can see from the comparison of the curves 
shown in Fig. [2] at the left. However, the value of the difference is relatively small, and, even for x = 2, it is on the 
order of seven percent. We also note that the strongest effect of nonlinearity, which makes the wave profile steeper, 
manifests itself as early as in the region of x where the difference between numerical and analytic results is small: for 
the example under consideration with ajv — 10, the corresponding value of the vx coordinate proves to be on the order 
of 0.08. To illustrate the possibility of increasing the accuracy of analytic calculations with the use of approximate 
symmetry, in Fig. [2] (right), we present the curves showing the difference between the zero-approximation solution 
q^ 1 and the solution q obtained from numerical simulation. One can see that the change from q(°) to qW already 
considerably improves the agreement between the numerical results and the approximate analytic ones. 



V. CONCLUSION 



Summarizing the results, we note that, in the present paper, we used the modified generalized Webster equation as 
the model to study exact and approximate analytic solutions to the problem of sound wave propagation in a channel 
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FIG. 1: Curves representing the solution to modified generalized Webster equation (|10[) versus quantity r. The curves are 
obtained from Eqs. (|33|) and (|35[1 for the periodic initial condition W(£) = cos£ for different values of the x coordinate along 
the axis of the channel whose cross-sectional area varies according to the exponential law \i(x)/v = exp(ai). The left-hand 
plot shows the curves for vx = 0, 0.2, 0.5, 1, 2, and 4 (the increase in vx corresponds to the passage from the upper curves to 
the lower ones at the axis r = 0) and for a/v = — O.f and a/v = f. The right-hand plot shows the curves for vx = 0, 0.08, 0.2, 
0.5, f , and 2 (the increase in vx corresponds to the passage from the upper curves to the lower ones at the axis r = 0) and for 
a/v — —0.1 and a/v = 10. 



with a varying cross section. The application of group methods allowed us to determine the specific types of cross 
section profiles for which the aforementioned problem allows new exact group invariant solutions. The application 
of renormalization-group symmetries allowed us to determine approximate analytic solutions for arbitrary initial 
conditions and sufficiently smooth profiles of variation in the waveguide cross section and to demonstrate the way of 
refining the approximate solutions. 

In addition to the conventional applications discussed in the Introduction, it is necessary to point out the possibility 
of using the results of our study in topical interdisciplinary investigations, which are primarily related to studies of 
biological tissues and medical applications. For example, equations of the type of the modified and generalized 
Webster equations had been used in hemodynamics for describing the nonlinear pulse waves fl8j ). Inhomogeneities 
of blood vessels are characterized by an inner cross-sectional area S(x) that decreases with distance from the heart. 
In addition, the stiffness of the vessel channel varies. The stiffness, which affects the wave propagation velocity, is 
determined by the varying ratio between the collagen and elastin contents. Peripheral vessels are usually stiffcr. The 
diagnostics of blood vessels is an important application of mathematical models based on the Webster-type equations 
and their complicated modifications, which are to be studied in the future. 
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of dimensional reduction and variable geometry effects in wave propagation problems described by equations of 
mathematical physics. 
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FIG. 2: Comparison between obtained solutions q to modified generalized Webster equation (|10|l by numerically solving Eq. 

(solid curves in the left-hand plot) and q = (solid 



(|10|l (dashed curves) and the approximate analytic solutions q 
curves in the right-hand plot) for different values of the coordinate: vx = 0.2, 0.5, 1, and 2. The increase in the ux coordinate 
along the pipe axis corresponds to the passage from the upper plots to the lower ones. The values of other parameters are the 
same as those in Fig. Q] 



